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ABSTRACT 

We consider the impact of quasar outflows on structure formation. Such outflows are potentially more important 
than galactic winds, which appear insufficient to produce the level of preheating inferred from X-ray observations 
of galaxy clusters. At late times, energetic material from the densest objects in the centers of galaxies makes 
its way into the intergalactic medium (IGM), impacting structures on many scales, much as supernovae impact 
structures on many scales within the interstellar medium. Using a simple analytical model for the distribution 
of quasars with redshift, coupled with a one-dimensional Sedov-Taylor model for outflows, we are able to make 
robust statements about these interactions. As large regions of the IGM are heated above a critical entropy of 
Sail ~ 100 keV cm^, cooling become impossible within them, regardless of changes in density. On quasar scales, 
this has the effect of inhibiting further formation, resulting in the observed fall-off in their number densities below 
z « 2. On galaxy scales, quasar feedback fixes the turn-over scale in the galaxy luminosity function {L-k) as the 
nonlinear scale at the redshift of strong feedback. The galaxy luminosity function then remains largely fixed after 
this epoch, consistent with recent observations and in contrast to the strong evolution predicted in more standard 
galaxy-formation models. Finally, strong quasar feedback explains why the intracluster medium is observed to 
have been pre-heated to entropy levels just above Smt, the minimum excess that would not have been erased by 
cooling. The presence of such outflows is completely consistent with the observed properties of the Lyman-a 
forest at z ~ 2, but is expected to have a substantial and detectable impact on Compton distortions observed in the 
microwave background and the multiphase properties of the "warm-hot" (z = 0) circumgalactic medium. 

Subject headings: quasars: general - galaxies: evolution - intergalactic medium - large-scale structure of the 
universe - cosmology: theory 



1. INTRODUCTION 

Non-gravitational heating is known to have played a central 
role in the formation of galaxy clusters. This is most clearly 
illustrated by the observed discrepany in the X-ray Luminos- 
ity Temperature (Lx - T) relationship of the diffuse intracluster 
medium (ICM). If non-gravitational heating was unimportant, 
the gas density distribution would be self-similar, resulting in 
Lx (X (Kaiser 1986), but instead the observed slope steep- 
ens considerably for low-temperature clusters (e.g. David et al. 
1993; Arnaud & Evrard 1999; Helsdon & Ponman 2000). Non- 
gravitational heating also manifests itself in the anomalous hot 
gas fractions observed in clusters [e.g. Muanwong et al. 2002; 
Tornatore et al. 2003), as well as the low overall level of the 
X-ray background (Pen 1999). While the origin of this heat- 
ing remains unknown, a wide variety of arguments point to the 
importance of active galactic nuclei, and in particular the lu- 
minous quasars. Although high-redshift starburst galaxies are 
observed to host massive supernova-driven outflows (e.g. Pet- 
tini et al. 2001), it is unlikely that they are able to generate 
the ~ 1 keV of energy per gas particle required to account for 
cluster preheating (Cavaliere, Menci, & Tozzi 1998; Balogh, 
Babul & Patton 1999; Kratsov & Yepes 2000; Wu, Fabian, 
& Nulsen 2000; Bialek, Evrard & Mohr 2001; Borgani et al. 
2001; Brighenti & Mathews 2001; Tozzi, Scharf, & Norman 
2001; Babul et al 2002). In fact, this level of heating is so high 
that it greatly exceeds the energy available from the Type-II 
supernovae resulting from metal-enriched stars, as well as the 
additional energy input from pair-production supernovae from 
a possible generation of massive primordial stars (e.g. Scanna- 



pieco, Schneider, & Ferrara 2003). Quasars, on the other hand, 
represent a vast reservoir of energy that is largely untapped in 
current theoretical modeling. 

Radiation from quasars should not be significantly coupled 
to the surrounding gas, except perhaps in the innermost cool- 
ing flow regions of the ICM where Compton heating or cool- 
ing could be important (Ciotti & Ostriker 1997, 2001). Nev- 
ertheless, the overall bolometric luminosities of quasars are far 
above the levels necessary to effectively preheat the intraclus- 
ter medium (e.g. Inoue & Sasaki 2001; Nath & Roychowdhury 
2002). Furthermore, a large fraction of quasars are observed to 
host massive outflows of material, with energies comparable to 
their observed luminosities. Thus the light observed to be ema- 
nating from quasars at intermediate redshifts, while not directly 
causing gas heating, is likely to be strongly correlated with out- 
flows that provide a substantial source of non-gravitational en- 
ergy. 

Like their high-redshift counterparts, these outflows would 
have consequences beyond their impact on the intergalactic gas. 
Detailed studies of high-redshift dwarfs indicate that starburst- 
driven shells are able to strip gas out of nearby overdense re- 
gions that would have otherwise later formed stars. This gas is 
ejected into space, leaving behind empty "dark halos" of non- 
baryonic matter (Scannapieco, Ferrara, & Broadhurst 2000). 
As it is much easier to strip gas from small objects, such winds 
naturally exert a feedback effect on neighboring dwarf galaxies 
(Scannapieco, Ferrara, & Madau 2002), causing a tremendous 
decrease in their numbers and slowing the formation of further 
outflows. 

The high entropy of preheated gas will also strongly affect 
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galaxy formation, as any level of entropy that can perturb gas 
density profiles in clusters will obviously have a major impact 
on the much shallower potential wells of galaxies. A moderate 
increase in entropy has been seen to reduce angular momen- 
tum loss during collapse in numerical simulations, (Mo & Mao 
2002; van den Bosch et al. 2003), and may explain the reduced 
X-ray luminosities of late-type galaxies (Mo & Mao 2002). 
Further increases in entropy are able to cut-off the supply of 
cold gas available for galaxy formation, resulting in strong neg- 
ative feedback, which may be important in determining galaxy 
morphology (Oh & Benson 2003), and which has direct impli- 
cations for clusters. 

The entropy to which clusters have been preheated is not 
arbitrary, but rather corresponds approximately to the thresh- 
old value for cooling within the age of the universe (Voit & 
Bryan 2001). Indeed, for S > Sent ~ 100 keV cm^, gas can 
never cool within a Hubble time regardless of the density to 
which it is compressed, and is therefore permanently unavail- 
able for galaxy formation (Oh & Benson 2003). Thus if non- 
gravitational heating in clusters had raised the gas to a some- 
what lower adiabat, this energy would have been quickly radi- 
ated away, and local X-ray observations would have left us none 
the wiser as to its presence. Yet, if clusters had been heated 
much above Sent, such observations would have easily detected 
this excess. While this coincidence can only result from fine- 
tuning in a model that accounts only for gas heating, it is an 
unavoidable feature of a model that accounts for its feedback 
on the process of galaxy formation itself. As quasar outflows 
heat the gas around them, they shut off gas cooling, both in the 
intracluster medium and in the gas that is accreting into smaller 
objects that would have otherwise formed into galaxies. In this 
case, the supply of cold gas available for star formation and 
quasar fueling is determined as a result of a strong feedback 
equilibrium reached between radiative cooling and quasar heat- 
ing. 

It is intriguing that the redshifts at which the cores of groups 
and clusters were assembled (z ~ 2 - 3) also correspond to the 
peak of quasar activity. Oh & Benson (2003) explored the ef- 
fects of global preheating at z ~ 2 - 3 on star-formation, cut- 
ting off the supply of cold gas, but did not consider the spatial 
dependence of entropy injection. In fact, more highly biased 
regions should be heated to high entropy at earlier epochs, im- 
plying that star-formation is self-limiting in the highest density 
peaks. A similar self-Umiting mechanism may also operate in 
high-redshift star formation (z > 10), when the residual entropy 
from early reionization prevents star formation in low mass ha- 
los (Oh & Haiman 2003). However, early preheating due to 
galactic winds at high redshiftz ~ 10 would have little effect on 
the present-day galaxy population, since the residual entropy of 
the IGM heated at high-redshift is too low (Benson & Madau 
2003). 

While several groups have recently studied the relationship 
between the formation history of galaxies and the cosmological 
evolution of quasars, none have accounted for the global impact 
of quasars on the IGM and its consequences for future galaxy 
formation. Kauffmann & Haehnelt (2000) incorporated a sim- 
ple scheme for the growth of supermassive black holes into a 
semi-analytic model of galaxy formation, reproducing the ob- 
served luminosity function of quasars in the context of major 
mergers. Ciotti & Van Albada (2001) used the observed rela- 
tionship between the masses of supermassive black holes and 
the velocity dispersion of their hosts to study the possibiUty of 
substantial dissipationless merging in the formation of ellipti- 



cal galaxies. Quasar formation regulated by local feedback was 
studied by Wyithe & Loeb (2003; hereafter WL03), as will be 
discussed in more detail below. Di Matteo et al. (2003) stud- 
ied the metalUcity evolution of quasar host galaxies, finding 
metallicity gradients and trends consistent with observations. 
Haiman, Ciotti, & Ostriker (2003) studied how the duty cycle 
of quasar activity can be constrained from the optical quasar lu- 
minosity function and the masses of host galaxies. Finally, and 
perhaps most directly related, Granato et al. (2003) studied the 
co-evolution of quasars and spheroidal galaxies, and the impact 
of quasar-driven winds on the gas within such galaxies. 

In this paper we adopt a more global view, in an attempt to 
understand the features that emerge when quasar feedback is 
accounted for in current scenarios of structure formation. Us- 
ing a simple analytical model for the distribution of quasars 
with redshift, coupled with a one-dimensional Sedov-Taylor 
model for the outflows associated with these sources, we ana- 
lyze the impact of quasars on galaxy formation, the structure of 
the intergalactic medium, and the further formation of quasars 
themselves. Despite the relative simplicity of our approach, we 
are able to make robust statements as to the general features of 
strong feedback models that differ from more standard scenar- 
ios, as well as indicate which observational constraints can be 
used to most easily discriminate between the two. 

Throughout this paper we restrict our attention to a Cold 
Dark Matter (CDM) cosmological model with parameters h = 
0.7, n„ = 0.3, Qa = 0.7, Qb = 0.045, = 0.87, and n = 1, where 
h is the Hubble constant in units of 100 km s"' Mpc"', fl^, n\, 
and are the total matter, vacuum, and baryonic densities in 
units of the critical density, cxg is the variance of linear fluctua- 
tions on the 8/z~'Mpc scale, and n is the "tilt" of the primordial 
power spectrum. These choices are motivated mainly by mea- 
surements of the cosmic microwave background, the number 
abundance of galaxy clusters, and high-redshift supernova dis- 
tance estimates (e.g. Spergel et al. 2003; Eke et al. 1996; Perl- 
mutter et al. 1999). Finally, we make use of the Eisenstein and 
Hu (1999) fit to the transfer function in computing the matter 
power-spectrum. 

The structure of this work is as follows. In §2 we describe a 
simple model of the quasar luminosity function (WL03), which 
relates the formation rate of quasars to the merger rate of col- 
lapsed dark-mater halos. In §3 we consider the scales of quasar 
outflows, and in §4 we develop a model for the spatial distri- 
bution of the resulting high-entropy material and its feedback 
on structure formation. In §5 we consider the impact of this 
feedback on the formation of quasars and galaxies, contrast it 
with higher-redshift feedback effects, and examine the impact 
of quasar outflows on the observed properties of the IGM at 
high and low redshifts. A discussion is given in §6. 

2. QUASAR POPULATION 

The first step in studying the impact of quasars on structure 
formation is to model their overall distribution. As discussed 
by Yu & Lu (2003), any such model is limited by two major 
constraints. The first of these is the observed number density 
of quasars, which has been well-measured out to redshifts ^ 6. 
Secondly, as quasars are believed to be powered by supermas- 
sive black holes (e.g. Blandford & Znajek 1977), such models 
must reproduce the measured black hole mass - bulge velocity 
dispersion (Mbh- Cc) relation. 

While a wide variety of models of quasar formation have 
been suggested within the CDM paradigm (Efstathiou & Rees 
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1988; Haehnelt & Rees 1993; Kauffmann & Haehnelt 2000; 
Nulsen & Fabian 2000; Monaco, Salucci, & Danese 2000; 
Haiman & Loeb 2001; Wyithe & Loeb 2002; Menci et al. 2003; 
Bromley, Somerville, & Fabian 2003), these can be roughly di- 
vided into two types. In early studies, quasar formation rates 
were tied to the formation rates of dark matter halos, as given 
by the rate of change in the Press & Schechter (1974, here- 
after PS) mass function of dark mater halos. As both observa- 
tional constraints and theoretical tools improved, such models 
were widely superseded by approaches that tied quasar forma- 
tion to mergers between dark mater halos. Theoretically, this 
had the advantage of eliminating poorly defined "halo forma- 
tion times," as well avoiding corrections due to the fact that 
the derivative of the PS mass function counts both the forma- 
tion of objects of a given mass as well as their "destruction" by 
merging into larger objects. Observationally, on the other hand, 
there is considerable evidence to support the idea that quasars 
are formed in mergers (e.g. Osterbrok 1993). A large fraction 
of quasar hosts display a disturbed morphology (e.g. Smith et 
al. 1986), have nearby companions (e.g. Vader et al. 1987), or 
even posses features that can be interpreted as tidal tails (Stock- 
ton & Ridgway 1991). Furthermore, the existence of the tight 
Mhh - Cc correlation suggests that quasar formation is closely 
related to bulge formation, which itself is thought to be trig- 
gered by mergers (e.g. Barnes & Hemquist 1992). Finally, we 
note that even the merger prescription is not enough to model 
the low-redshift quasar population, which exhibits a strong de- 
crease at z <J 2. This is traditionally interpreted as due to the 
ever-increasing gas cooling times within larger halos (e.g. Rees 
& Ostriker 1977), but also arises in a model of strong quasar 
feedback, as we shall see below. 

As the focus of our investigation is not the detailed model- 
ing of the processes of quasar formation, but rather the global 
feedback of this population on structure formation as a whole, 
we do not attempt a survey of these quasar models. Rather we 
make use of the simple and elegant model developed in WL03 
and Wyithe & Loeb (2002), which satisfies both types of con- 
straints in the context of major mergers. Furthermore, as the 
model is constructed only to reproduce the distribution of z ^ 2 
quasars, it contains no prescription for gas cooUng and accre- 
tion at low redshifts that might otherwise obscure our results. 

Instead, the authors simply assume that following a merger, 
a black hole shines at its Eddington luminosity (as observed by 
e.g. Willott, McLure, & Jarvis 2003), and returns a fraction of 
this energy to the galactic gas, eventually disrupting its own fuel 
source. In this case a quasar shines for the dynamical time of the 
cold gas surrounding the black hole, which is usually located in 
a disk with a characteristic radius (A/ \/2)rvir ^ 0.035rvir (Mo, 
Mao, & White 1998). The corresponding dynamical time is 
then (e.g. Barkana & Loeb 2001) 



tayn ~ 0.035rvir/vc w 5 x lO'yr (1 +z)' 
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(1) 



where A ^ 0.05 is the spin parameter, rvir is the virial radius, 
and Vc is the circular velocity of the galactic halo, which can be 
written as 



Vc = 140 km s"^M!f(l+z) 



^1/2 



(2) 



where Mi 2 



Afhalo 



Note that the timescale given by eq. 



(1) is comparable to the e-folding time of the black hole mass 
or Salpeter time ts ^ 2.7 x 10^(erad/0.06)yr, where Crad is the 
radiative efficiency. 



Furthermore, this "self -regulation" condition amounts essen- 
tially to assuming that the total energy output from a quasar is 
proportional to the depth of the potential well of the halo in 
which it is contained, implying Mbh oc v^. Note however that 
this index is slightly larger than the 4-4.5 inferred from the lo- 
calMbh-cr^ relation (Merritt & Ferrarese 2001; Tremaine et al. 
2002), the difference having its origin in the observation that 
the Vc - relation is shallower than linear (Ferrarese 2002). In 
this model the Mbh - Vc relation is independent of redshift as 
observed by Sheilds et al. (2003), with 

^- = ^-4^1«^^«^(300i^)'- 
Combining this relation with eq. (2) gives 

Mbh(Mhalo , Z) = ebh(Mhalo , z)Mhalo , (4) 

where ebh(Mhdo,z) = ebh.oMj^ 1 +z)2 , and ebh.o = 3.2 x 10""^^. 
The only adjustable parameter in this model is F, which, in 
our notation, was taken to be 0.6 in WL03, and is taken 
to be 1.0 in our fiducial modeling below. Observations 
(e.g. Ferrarase 2002) suggest that Mbh = M0(1.66 ± 0.32) x 

10^ [cr,/(200kms-')]'*-^^^°-^^ with logio(vc/300) = (0.84 ± 
0.09) logio(o-c/200km/s) + (0.0 ± 0.19), which would aUow for 
models with F values between 0.15 and 7.0 to fall within la- 
experimental errors, due primarily to the uncertainty in the re- 
lation between and (Tc- 

2. 1 . Optical Quasar Luminosity Function 

Having developed a model for quasar formation during 
mergers, we are now able to construct a luminosity func- 
tion, '^(z,Lb), which measures the comoving number density 
of quasars per unit B-band luminosity (Lb). Taking the ob- 
served relation of Lboi = IQ.ALb (Elvis 1994) means that Mbh = 
^Edd,B/(5.73 X IO^LqMq), and we can invert eq. (4) to com- 
pute Mhaio as a function Lb, relating ^(z,M) to the merger rate 
as 



fdyn(z) 



5ebh 5.73 X lO^LgMr' 



O.SMhalo 



C^AMhalo (5) 



X^fc^'):^5^fcAMha,o,Z,Mhalo), 



<iAMhalo<i? 



where M' = Mhaio- AMhaio- Here the comoving density of halos 
with mass M is given by the usual PS prescription 



dn p 



'dM 



(M) 



f[u(z),a\M)'\, 



(6) 



where p is the mean matter density, ^(z) = l.69D~^(z) is the 
"collapse threshold" at a redshift z with D(z) the linear growth 
factor, cr^(M) is the the variance of linear fluctuation within 
spheres containing a given mass M, and 



f(iyy)da^ 



exp 



2a2 



da^ 



(7) 



is the mass fraction contained in halos with masses associated 
with variances in the range from to a^+da^. Finally 

(Zp,Mp,Zf,Mf)=-^^—f2 [up,a\Mp),Vf,a^(Mf)\ 



dMpdt 



dt(zf) dvf 



(8) 
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is the rate of mergers of halos with masses between Mp and 
Mp + dMp and formation redshifts Zp into halos of mass Mf at a 
merger redshift Zf and 



4^' 



f{Vp,(Tl) 



(9) 



is the conditional probability that an object with a mass Mp at 
a redshift Zp will be found in an object with a mass My at zj 
(Lacey & Cole 1993). 

The resulting luminosity functions are shown in Figure 1. 
Again we note that the WL03 model is constructed to match the 
quasar distribution at high redshifts, and thus assumes that all 
the gas within a galactic halo cools on a time much shorter than 
the Hubble time. This is the primary cause of the overestimate 
of at the lower redshifts, as eq. (6) includes an additional 
contribution from massive halos, which form into groups and 
clusters rather than individual galaxies. 

3. QUASAR OUTFLOWS 

In order to relate '^(z- Lb) to the overall level of kinetic feed- 
back, we must adopt a simple model for the formation and prop- 
agation of quasar outflows. Here two major types of objects 
contribute. In radio-loud (RL) quasars, the outflow takes the 
form of a coUimated jet, which deposits particles into a cocoon 
that expands into the surrounding medium. These jets, which 
are found in ~ 10% of all quasars, are observed to have over- 
all kinetic luminosities {Li) that are correlated with the bolo- 
metric luminosity of the quasars, with 0.05 ^ Lk/Lboi ^1-0 
(Willott et al. 1999). The major radio-quiet (RQ) sources, on 
the other hand, are broad absorption line (BAL) quasars, whose 
optical absorption troughs are thought to correspond to outflow- 
ing clouds with velocities up to O.lr. Similar velocities are 
also observed in outflowing material detected at X-ray wave- 
lengths {e.g. Pounds et al. 2003). Although these make up ap- 
proximately 10% of the observed population of quasars, it is 
believed that all quasars host these outflows, with a 10% cov- 
ering factor {e.g. Weymann 1997). Again observational and 
theoretical arguments suggest that the impact of these clouds 
on the surrounding medium can be parameterized with overall 
kinetic luminosities correlated with the bolometric luminosities 
with 0.\<,Lk /Ltoi <, 1 (Furlanetto & Loeb 200 1 ; Nath & Roy- 
chowdhury 2002). 

Do quasars that host outflows differ widely in their proper- 
ties from the general quasar population? In the radio-loud case, 
this relation is a subject of debate. A number of earlier stud- 
ies claimed that, with substantial overlap, RL quasars tended 
to be found in Abell 0/1 clusters, while RQ quasars generally 
inhabit less-dense regions {e.g. Yee & Green 1984; EUingson 
1991). More recent Hubble Space Telescope studies, however, 
have shown that to some degree this is a selection effect, with a 
RQ population being found in clusters {e.g. Bachall et al. 1997), 
perhaps in the same relative proportions as RL quasars (McLure 
& Dunlop 2001). 

In the BAL case, on the other hand, the primary issue is not 
if the spatial distribution of objects is different from the general 
quasar population, but rather whether the presence of broad ab- 
sorption lines represents a stage in the life cycle of quasars {e.g. 
Briggs et al. 1984), or is simply a geometric effect {e.g. Wey- 
mann et al. 1991). In this case the geometric interpretation is 
favored by the similarity between the optical continuum and 
Une emission of BAL and non-BAL quasars (Boroson & Mey- 
ers 1992), as well as by the fact that the sub-millimeter emission 



of quasar host galaxies is not related to the presence of broad 
absorption lines (Willott et al. 2003). 

Finally there is the issue of outflow structure, which is clearly 
jet-like in RL quasars, but may later expand roughly spherically 
(Furlanetto & Loeb 2001). In BAL quasars, outflows are likely 
to instead take the form of axisymmetric winds, accelerated by 
an unknown mechanism such as magnetrocentriful forces {e.g. 
Blandford & Payne 1982) or radiation pressure from lines {e.g. 
Proga, Stone, & Kallman 2000). 

For a review of the observed relations between RL, BAL, and 
other quasars, the reader is referred to Antonucci (1991). For 
the purposes of this study however, the important point is orily 
that the populations of quasars that are observed to host out- 
flows are consistent with a random subset of all quasars, and 
that these outflows exhibit a variety of structures. In this case, 
we adopt the simplifying assumption that all quasars host out- 
flows of material, with a kinetic energy input equal to some 
fraction, Ck ~ 0.05, of the total bolometric energy released over 
the quasar lifetime (see §5.2 for justification of our choice of 

Furthermore, the source hfetime fdyn ~ lO^yr is generally 
shorter than other time scales under consideration. We can then 
approximate the energy injection as a point source explosion of 
energy £ k, ejjihoifdyn- We shall assume the bubble evolves adia- 
batically and ignore energy losses due to radiative cooling, pdV 
work against the IGM, and gravity. We also ignore the Hubble 
expansion and peculiar velocities (including accretion infall), 
and assume that the bubble expands into a medium of uniform 
overdensity 5s. The effect of relaxing these assumptions will be 
exanained later. 

3.1. Ou^ow Scales 

Under these simpUfications, the bubble behaves as a Sedov- 
Taylor blast (Sedov 1959) solution, which describes the adia- 
batic expansion of a hot sphere of plasma into cold medium. 
For our cosmological model, this gives a physical radius of 

/e;,2x1/5 

= 6 (^— j =1 -VMpc^/' 5-'" (1 +z)-3/5 (10) 

where E(,q is the energy in the hot medium in units of 10^° 
ergs, feyr is the expansion time of the bubble in Gigayears, 
the overdensity of the surrounding medium 5s is defined such 
that 5s = Pb/Pb (rather than Pb/ph- 1) and = 1.17 for a 
7 = 5/3 gas {e.g. Shu 1992). Thus the shock velocity is 

3/2 E^l/2.-l/2 



2/5R/t = 1500/?:-;£,„ 



5s (1+z) "^^^ km s where 



^s,Mpc = ^sMpc ^ and, assuming a strong adiabatic shock, the 
postshock temperature is 7^ = {3fimpV^)/{l6kB) = 13.6^rv^/(km 



s ) . We therefore obtain the temperature as a function of 
proper radius 



T; = 3.1 X lO'^KEeod:^ (1 +z)"^/?r^ 



SjMpc • 



(11) 



The functional form of this equation can be easily understood 
from E oc pRl{kBTs)= constant, i.e. the blast wave conserves 
energy. As the postshock density psh = {l+^)/{j- 1)/5igm<5s = 
4yO/GM<5s, for 7 = 5/3, the postshock entropy can then be calcu- 
lated as 

Ss{r) = 4« = 1-8 X lO^keVcm^Efio-^:'/' (1 +z)"'/?;,Mpc- (12) 



2/3 
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Finally can be related to the halo mass and merger redshift 
corresponding to a given quasar, by combining eqs. (1) and (4), 
yielding 

(13) 



E6o = l.2ei,FM%\l + z) 



4. CRITICAL ENTROPY 

The simple model of entropy injection described above can 
be easily extended to capture the global spatial dependence of 
energy injection and coohng, which has not been considered 
in previous preheating models. The great advantage of con- 
sidering entropy is that it is conserved during adiabatic pro- 
cesses such as subsonic compression and the Hubble expan- 
sion, and remains constant as long as radiative cooling is in- 
efficient. Moreover, Oh & Benson (2003) pointed out that for 
5(1 +z)^ > 10, the cooUng time depends only weakly on den- 
sity and is always greater than the Hubble time. Thus, gas that 
is shocked to entropies greater than a critical value can never 
cool, regardless of the densities to which it is compressed. 

This can be most clearly understood if we write the cooUng 
time as a function of temperature and entropy, rather than tem- 
perature and density. Using S = T/n^^^, the isobaric cooling 
time can be written as 



^cool ■ 



i3/2)nkBT 



3 //i^ 

2\ii ) T^I^K{T) 



(14) 



where for a fully ionized gas, fi = 0.62, and /i^, = 1.18. The 
quantity in square brackets is a function of temperature only, 
and the entropy only changes the overall normalization of the 
coohng time. 

In the left panel of Figure 2, we plot the cooling time at a 
constant entropy as a function of temperature for Z = 0. IZ© and 
O.3Z0, and 5ioo = 5/(100 keVcm^) = 1. We use fits to the cool- 
ing function computed by Sutherland & Dopita (1993). Chang- 
ing 5ioo shdes the curves up and down the vertical axis, but 
does not change their shape. The cooling time has a deep min- 
imum at Tmin = 2.3 X lO^K, where the cooling function peaks. 
Although this time can be shorter than fcooKTmin) at high tem- 
peratures, when A(r) oc r'/^ due to free-free emission and 
^cooi oc this only occurs if T > lO^K. As such hot gas can 
only be retained in the potential wells of extremely massive 
clusters with Tyir > 10 keV, we can safely ignore this regime. 

Thus, if fcooi(7min) > tn, the gas can never cool in a Hubble 
time, regardless of any subsequent expansion or compression. 
Since ?cooi(7min) oc 5^/^, this will occur for all entropies greater 
than some critical value San, given by 



5crit = 280keVcm2(l+z)-i 

A(7;„in) 



(1+2)3/2 



2/3 



6.3 X 10~2^ergs~icm3 



2/3 



(15) 



where E{z) = ^fim(l+z)'^ + ^^A- We plot this critical entropy, 
expressed in terms of -Sioo.crit = 5crit/(100keVcm-^) in the right 
panel of Figure 2. 

The fact that the entropy floor seen in groups and clusters is 
greater than 5crit is no coincidence: all lower-entropy gas would 
have dropped out by radiative cooling. Indeed, it was this very 
fact which lead Voit & Bryan (2001) to suggest that radiative 
coohng could cause the observed entropy floor, although cool- 
ing alone could not be wholly responsible, since in the absence 



of preheating this would require massive overcooling in groups 
(e.g. Balogh et al. 2001; Oh & Benson 2003), and would result 
in entropies at large cluster radii that are much lower than ob- 
served (Voit & Ponman 2003). The fact that, at all radii, the 
observed entropy floor is not far above 5crit, however, is likely 
to be a manifestation of strong feedback, as discussed in detail 
below. 

4.1. Spatial Distribution 

The presence of a critical entropy suggests that the IGM can 
be divided into two fractions, material that has been heated to 
above 5crit and is therefore unavailable for galaxy formation, 
and a cold phase at T ^ 10'*K where all memory of preheating 
has been erased. At this point an analogy with the well-known 
(McKee & Ostriker 1977) three-phase model of the intersteUar 
medium (ISM) within galaxies is irresistible. In the ISM case, 
supernova explosions result in a three-component medium in 
which a large fraction of the Galactic volume is filled with hot, 
tenuous gas. The remainder of the cold gas is then divided into 
dense clouds in which star formation occurs, and a relatively 
cool intercloud medium that accretes onto the clouds that lie 
outside of hot bubbles. 

In the cosmological case, the 5 > 5crit phase, the T ~ 10'' 
phase, and the galaxies themselves play the roles of the super- 
nova bubbles, the cool medium, and the dense clouds, respec- 
tively. Quasar outflows expand into the IGM, leaving behind 
heated regions in which no further structure formation can oc- 
cur. In the cooler zones between these bubbles, gas accretion 
continues apace, resulting in the formation of continuaUy larger 
structures. Finally, the gas within the galaxies themselves not 
only condenses to form stars, but also fuels supermassive black 
holes, seeding new bubbles of 5 > 5crit intergalactic material. 

We are interested in the radius /?heat of the bubble that is 
heated above the critical entropy, S > Sent. To find this, we con- 
sider the sequence of adiabats through which the gas is shocked 
and invert eq. (12). This gives a proper radius of 

/?heat = 5.6Mpc5ioo,crit(z)"'/'£6o''5s"''"' (1 +zT'^' , (16) 

(17) 

and a corresponding gas mass Mb.heat = (47r/3)/?^g^(/9/,5s of 

Mb,heat(^,Z,M) = 4.6 X lO^^Mg 5ioo,cri,(z)"' £60 ^s"'^^ (1 +2)"' • 

(18) 

If we equate this to the total shocked mass as a function of time, 
Msit), via equation (10), we can find the timescale on which 
this entire mass of gas is preheated. It can be quite long for 
low-redshift, energetic sources: 

^expand = 20Gyr5ioo,crit(z)"'/%'o^'^"'/'(l +z)-^'\ (19) 

Thus, even after the demise of the quasar population at low 
redshift, their hot bubbles continue to expand and encompass 
larger tracts of the IGM. We discuss this further in section 5.5. 

Unlike the ISM case, the initial positions of these sources 
are not determined by baryonic physics or galactic dynamics, 
but rather by the underlying dark-matter distribution. In order 
to track the evolution of the hot baryonic component, then, we 
must combine our outflow modeling with the formation rate of 
massive dark matter halos. Thus, the average number of sheUs 
with S > Sent impacting a random region of space with an over- 
density of 1 , at a redshift z can be computed by integrating over 



6 



the range of masses that contribute to the overall luminosity 
function at each redshift. This gives 



iz)} = — dz—i 
Pbfi Jz dz' 



dyn 



Mb,„j„(z', l,z,M) -^^[z' ,Lb{7! ,M)'\, (20) 
dM 



where p^.o is the mean baryonic density at z = 0, M4o(z') is the 
mass corresponding to a circular velocity of 40 km s"' , which 
approximates the minimum mass that can form after reioniza- 
tion {e.g. Barkana & Loeb 2001), and 



M, 



b.min 



Mb.heat(l,z',M), —p^5sRl(z,z,M) 

(21) 

accounts for the possibiUty that a given outflow may not have 
expanded to contain a mass Mb.heat by the final redshift z- The 
factor of fjyj, accounts for the fact that the finite duty cycle of 
quasars implies that only a fraction ~ (fdyn/fn) are visible at any 
given time; the overaU number density of sources is a factor 
~ (fH/fdyn) higher than the observed number density of bright 
sources. Assuming that outflows are randomly distributed al- 
lows us to simply relate {Nsiz)) to the total mass fraction of gas 
shocked above 



nsiz) = l-sxp[-{Ns(z)}]. 



(22) 



As our model associates a mass and formation redshift with 
each quasar, it can be easily generalized to account for the in- 
creased number of 5 > 5crit outflows that impact the regions that 
later form into large galaxies and clusters, which is due to the 
fact that quasars tend to be biased towards the densest regions 
of space. In this case, the average number of such outflows im- 
pacting a region that forms a halo of mass Mf at a redshift Zf 
is 



{Nsizf,Mf)) = 



1 



d^ 



,dt\ 
dz' 



.1 Myn 



dMx 



pbfi Jzf ■ JMm(z') 

Mb 

,min 



(23) 



b{z' ,M,Zf,Mf)-^nz! Xb{z! .M)\. 
dM 

where is a now a "bias factor" that accounts for quasar clus- 
tering and Ssiz',Zf) is the overdensity of the forming object at 
the time it encounters an outflow. 

To compute b we employ an analytical formalism that tracks 
the correlated formation of dark matter halos. In this model, 
described in detail in Scannapieco & Barkana (2002), objects 
are associated with peaks in the smoothed Unear density field, 
in the same manner as the standard PS approach. This for- 
malism extends the standard method, however, using a simple 
approximation to construct the bivariate mass function of two 
perturbations of arbitrary mass and collapse redshift, initiaUy 
separated by a fixed comoving distance (see also Porciani et al. 
1998). From this function we can construct the number density 
of source halos of mass M that form at a redshift z at a comov- 
ing distance r from a recipient halo of mass Mf and formation 
redshift z/: 



dn , d^n 

^(M,z,r\Mf,Zf)=^^^iM,z,Mf,Zf,r) 



dn 
dMf 



iMf,Zf) 



(24) 



where -^(MfjZf) is the usual PS mass function and 

d^Mj(.M,z,Mf,Zf,r) is the bivariate mass function that gives 

the product of the differential number densities at two points 
separated by an initial comoving distance r, at any two masses 
and redshifts. Note that this expression interpolates smoothly 
between all standard analytical limits: reducing, for example, 
to the standard halo bias expression described by Mo & White 
(1996) in the limit of equal-mass halos at the same redshift, and 
reproducing the Lacey & Cole (1993) progenitor distribution in 
the limit of different-mass halos at the same position at differ- 
ent redshifts. Note also that in adopting this definition we are 
effectively working in Lagrangian space, such that r is the ini- 
tial comoving distance between the perturbations. Fortunately 
as shock propagation is likely to be more closely dependent on 
the column depth of material between the source and the recip- 
ient than on their physical separation, this natural coordinate 
system is more appropriate for this problem than the usual Eu- 
lerian one. We then define the bias factor as 



dn 

b{z,M,Zf,Mf) = ^(z,M,r\zf,Mf) 



(25) 



which is simply the number density of mass M, redshift z ha- 
los near the final object, divided by the average density of such 
halos. Finally for the effective radius, one works in Lagrangian 
comoving coordinates, such that r is the virial radius plus the 
radius which would contain Mb^min in a mean-density medium. 



47rpi,.o 



Mb.n 



1/3 



+ rvir(M). Note that this calculation of 



bias, while carried out within our merger formahsm, depends 
only on the mass of the quasar host galaxy, the turn-on red- 
shift for the quasar, and the mass and formation redshift of the 
final object. As these are the same quantities that arise in an 
approach based instead on accretion, one would expect the re- 
lationship between (A(s(z/,M/)) and (NsiM)) to be similar in 
any such model. This means that the rough features derived be- 
low are likely to persist, regardless of the detailed physics of 
quasar generation. 

To approximate (5s, we appeal to the standard top-hat collapse 
model (see e.g. Peebles 1980), which relates the true overden- 
sity to the "linear" overdensity, Sl expressing both quantities 
parametricaUy in terms of a collapse parameter, 0, as 



and 



Ss(z,Zf) = 



SLiz,Zf)=l + l ' (^-sin^)2/3 = 1.69^, (27) 

where Z)(z) is again the linear growth factor. 

Finafly, we relate {N(zf,Mf)) to fs(xf,Mf), the fraction of 
objects with collapse redshifts z/ and mass Mf whose formation 
is suppressed by preheating by quasar outflows. As our formal- 
ism can only account for the correlations between two halos, 
we must adopt an approximation for the fraction of objects that 
are affected by multiple outflows. The simplest approach is to 
assume that while outflows only impact a fraction of the overall 
cosmological volume, their arrangement within that fraction is 
completely uncorrelated (as in Scannapieco, Schneider, & Fer- 
rara 2003). In this approximation /5(z/,M/) and {Nsizf,Mf)) 



9 (6- sine f 
2(1 -cos 61)3' 



(26) 
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share the same relation as ilsiz) and {Ns{z)) in the average case, 
namely 

/s(z/,M/)= l-exp[-(Afs(z/,M/))]. (28) 
We will adopt this Ansatz throughout this paper. 

4.2. Theoretical Uncertainties 

Before we apply the formalism developed above to derive the 
features expected in a model of strong feedback, we first pause 
to assess the uncertainties involved in our approach. Primary 
amongst these is our use of an adiabatic Sedov-Taylor solution 
for a point source explosion in a uniform medium. Here we 
note that 

W « ^ = 4.8 X IQ^'^yrSl^^^E'isfl' (1 ^z)-'" , (29) 

Vsh 

which is larger than our source lifetime of ~ lO^yr (as in 
equation 1) justifying our assumption of instantaneous energy 
injection. 

Secondly, we ignored radiative losses. Since we are only 
interested in the early stages of the blast wave, when gas is 
shocked to high entropies (S > Send and can never cool, this is 
justified; /?heat is by definition the radius at which the cooling 
time of the gas is equal to the Hubble time. Beyond /?heat> ra- 
diative cooUng indeed becomes important and the evolution of 
the blast wave is no longer adiabatic. 

Note that Compton cooling will dominate over colUsional 
cooling at high redshifts (z > 7); however, we do not consider 
such high redshifts in our calculations: both the high IGM den- 
sities and the effectiveness of Compton cooling mean that the 
residual entropy is low. We also ignored the deceleration of the 
shock due to the pdV work it does on the IGM, which is jus- 
tified as long as pv^^ » ^gm, or Tshock > Tigm ~ IC'K, which 
is certainly the case. Gravitational deceleration is unimportant 
as long as the work done in sweeping up the shell out of the 
potential well W = 0.5 J pg^d^x = GTrMhaio^PiGM^sh is smaller 
than the energy of the explosion, which it is 

"^^^^ =4.5 X 10-5XA^/^^-/^1+z)-/^M,. 

(30) 

Note that since E oc Mj2^, W/E cc mJj*", so the dependence on 
halo mass is very weak. Finally, ignoring Hubble expansion is 
valid as long as Vsh H{z)Rsh, which is only marginally valid: 

^(Z)nheat 

Our model is therefore conservative in that we underestimate 
the extra kinetic energy associated with the Hubble expansion. 

In order to use the Sedov-Taylor solution, we have assumed 
a point source explosion in a homogeneous medium, character- 
ized by some mean overdensity S. In fact, the density profile is 
more likely to be that of a power law; in addition the medium 
is not static but is being accreted by the halo. It might seem 
therefore that the dynamics of the explosion are no longer ac- 
curately described. In fact, Barkana & Loeb (2001) and Furlan- 
etto & Loeb (2001) have constructed numerical models of out- 
flows and solved the coupled system of differential equations to 
describe the outflow dynamics, modeling both the radial depen- 
dence of the density profile as well as accretion infall. They find 
remarkable agreement with a simple Sedov-Taylor type model 
such as that used here, which ignores the density profile and 



infall, as well as cooling, external pressure of the IGM, and the 
gravitational potential of the host halo (see Fig. 1 of Furlanetto 
& Loeb (2001), and associated discussion). This can be under- 
stood, since in all cases /?heat ^ ^vir, and the swept-up mass is 
dominated by that accumulated at large radii, when the density 
is close to the cosmological mean density. It is therefore largely 
independent of the details of the infall region. 

We have completely bypassed the issue of how the quasar 
wind escapes from the host galaxy. Such issues are beyond 
the province of our semi-analytic model and can only be stud- 
ied with hydrodynamic simulations. To our knowledge, there 
has been no detailed numerical study of quasar winds, although 
there is a significant body of literature on starburst winds (e.g., 
Strickland & Stevens 1999; MacLow & Ferrara 1999; Mori, 
Ferrara & Madau 2002; Fujita et al. 2003). Simulations of a 
centrally concentrated nuclear starburst most closely approxi- 
mate quasar winds and generically find that starbursts punch a 
hole through the ISM of the galaxy ("blow-out"), venting most 
of the energy of the wind (with little radiative losses) while 
entraining little cold gas, which remains in the galaxy. Some 
fraction (few x 10%) of the kinetic energy of the outflow can 
be radiated in internal, inward propagating shocks if there are 
significant off-nuclear supernova explosions (Mori, Ferrara & 
Madau 2002), but the radiated energy diminishes for a single, 
central energy source. 

In fact, a single source, rather than smoothly distributed en- 
ergy injection, is the most efficient means of quickly producing 
blow-away (Strickland & Stevens 1999). The coupling between 
the ISM and the wind diminishes rapidly for more energetic 
bursts, as blow-out becomes more efficient; only for very small 
dwarfs Mg < 10^ Mq is a large fraction of the ISM entrained in 
the outflow (MacLow & Ferrara 1999). Blow-out would likely 
be particularly efficient for an initially bi-polar outflow such 
as an quasar wind. Thus, while an quasar wind has not been 
directly simulated, it certainly lies in the region of parameter 
space where energy should be efficiently delivered to the IGM. 
In any case, our model is only applicable once the wind has 
propagated out of the host galaxy: all uncertainty in energy loss 
within the host galaxy is encapsulated in the free parameter et. 
For simplicity, we have assumed Ck to be independent of galaxy 
mass. We emphasize that our model only taps a very small frac- 
tion, ek ~ 5%, of the total available kinetic luminosity, and so 
preheating can still be important even if an appreciable fraction 
of energy is radiated away in the ISM of the host. 

We approximate the metallicity of the IGM as O.IZ©, con- 
sistent with the level of pre-enrichment necessary to explain 
the lack of low-metallicity G-dwarf stars in the solar neigh- 
borhood {e.g. Ostriker and Thuan 1975), as well as in other 
large nearby galaxies (Thomas, Greggio, & Bender 1999). This 
should also be representative of the metallicity of regions sur- 
rounding quasar outflows at z = 2, when preheating by outflows 
peaks. Such regions now constitute the intracluster medium, 
which are observed to have metallicity Z w O.SZq; if 30% 
of all present-day stars have formed by z ~ 2, as indicated by 
the Madau plot, then Z(z = 2) ~ O.IZ©. However, assuming a 
higher metallicity does not greatly alter 5crit (see Fig 2). 

The greatest uncertainty is = L^/Lbou the fraction of the 
total bolometric luminosity (assumed to be the Eddington lumi- 
nosity) which appears as kinetic luminosity. Most observations 
suggest that Lk ~ Lb for both radio-loud (Willott et al 1999) 
and BAL quasars (Furlanetto & Loeb 2001), albeit with at least 
an order of magnitude uncertainty. Since ~ 0. ILboi (Elvis et 
al 1994), Ck ^ 0.1 is reasonably motivated from observations. 
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However, the uncertainty is so large that we shall simply take 
£/. to be a free parameter, which we use to match the observed 
quasar luminosity function, as described in the following sec- 
tion. 

5. IMPLICATIONS FOR STRUCTURE FORMATION 
5.1. Global Features 

From eqs. (24) and (28), we can construct the fraction of 
galaxies whose formation is suppressed given an arbitrary lu- 
minosity function ^[z,Lb(z,M)]. Yet our model is not com- 
plete without correcting ^ for the corresponding suppression 
in quasar formation. Naively, one might consider convolving 
the luminosity function with the fraction of objects affected 
by heating at a given redshift, replacing \E'[z,Lb(z,M)] with 
[1-/s(z,^] X ^{z,L^{z^M)]. Yet this approximation would re- 
sult in a severe underestimate, as in reality two galaxies of mass 
M/2 that formed at early times can still merge to ignite a quasar 
at z, even if the formation of additional galaxies of mass M is 
largely excluded at this redshift. 

A better approach relies on the fact, as we shall see be- 
low, that {Ns(z,M)) is in general monotonically increasing with 
time. Thus we define an "exclusion redshift," ZexiM) such that 
all M mass galaxies forming by ZeniM) contribute to the quasar 
luminosity function, while later forming galaxies do not, with 
Zex defined implicitly as 



(A^sUex(M),M]) = l. 



(32) 



With this definition, we then impose that for all halo masses and 
redshifts such that z < Zex(^^haio/2), eq. (6), gas can no longer 
cool and the supply of cold gas is shut off. Thus, only halos 
that formed at z > Zex will contain cold gas to fuel quasars, and 
only their mergers need be considered. These progenitor galax- 
ies will in general reside in a much larger host halo (such as a 
group or cluster), whose mass greatly exceeds the sum of the 
merging galaxy masses. Thus eq. (6) must be rewritten to in- 
clude an additional integral over possible final masses, M^ig of 
the final object. Defining M^jg = Mbig - AMhaio this integral can 
be written as 



5ebh5.73 X IO^LqM^ 



O.SMhalo 



(/ AMhaio X 



dn 
dM' 



-jiZe: 



,M') / dMl, 

JM' 



0.25Mhalo 



i^iz,.),cj\M')Mz),cT^i, 



d^Nr 



merge 



dAMbaiodt 



(Zex, AMhaio, Z, Mbig). 



(33) 



Note that here we assume that the galaxies contained in each 
sub-halo merge instantaneously as they fall into the larger halo 
Mhaio, and do not attempt to account for dynamical friction 
or tidal stripping within this region (e.g. Kauffmann, White, 
& Guiderdoni 1993; Cole et al. 1994; Somerville & Primack 
1999) as these are beyond the scope of our simple model and 
would only obscure the processes in which we are interested. 
Note also that in the limit in which Zc — > z, fi approaches 
(5[f7^(M')- cr^i'g] and thus eq. (33) reduces to eq. (6). 

We then adopt an iterative approach to obtain "^(z^Lb) and 
{Nsiz,M)) as a function of our feedback parameter e^: first 
computing (A(s(z,M)) from the simple WL03 model, apply- 
ing eq. (33) to correct the luminosity function, computing 



{Ns(z,M)) from this corrected function, and so on. In prac- 
tice this convergence is very quick, arriving at a self-consistent 
solution within two cycles, although we iterate five times in the 
results presented below. 

In Figure 3 we show the {Ns(zf,Mf)) values that result from 
models with F = I and three different values of e^. Here we 
choose Ek = 0.05 as our fiducial model, simply because it pro- 
vides the best fit to the quasar luminosity function, as discussed 
in §5.2. We then double and halve this value to examine the 
impact of varying the feedback efficiency on our modeling, re- 
sulting in the ek = 0.025 and Ck = 0.10 models depicted in the 
lower panels. In all cases {Nsizf,Mf)) is a monotonically de- 
creasing function of redshift, as assumed. 

While in a self-similar, hierarchical model, preheating by 
progenitors would be independent of the final host halo 
mass scale, three features of quasar outflows break this self- 
similarity. As energy injection is a non-linear function of halo 

5/3 

mass, E oc Mbh cx Mj^^^j^, feedback is a strongly increasing func- 
tion of mass. Secondly, the presence of a critical entropy im- 
poses a scale into the problem, as the number of supercritical 
outflows is self limiting, and thus does not greatly exceed one 
even for the largest objects. Finally, high mass objects are more 
strongly clustered, and therefore more likely to collapse out of 
preheated gas than low mass objects. The effect is similar to 
well-known fact that gas metallicity has a strong density depen- 
dence, and metals are likely to be clustered around the highest 
density peaks. Metals are likely to be a good tracer of gas en- 
tropy, and the two should exhibit similar spatial bias. 

Thus for larger masses halos with My ^ lO'^M©, {N(zf,Mf)) 
is a relatively weak function of mass. This is because as 
{N(zf,Mf)) > 1, the further formation of galaxies is strongly 
suppressed, and these objects are heated primarily by high- 
redshift M <C Mf progenitors. In this limit the bias is well ap- 
proximated by the Lacey & Cole (1993) progenitor distribution 



biM,z,Mf,Zf) = 



f[iyiz)-Hzf),aHM)-a\Mf)] 
f[p(zWiM)] 



(34) 



where i'(z) and i^(z/) are the "collapse thresholds" and a^(M) 
and a^iMf) are the variances associated with the progenitor 
and final objects, respectively. As M ^ My, a^(M) » cr^(M/), 
and therefore eq. (34) is nearly independent of the Mf. In 
other words, the number density of M ^ 10 '^Mq progenitors 
is roughly equivalent for groups and massive clusters that form 
at the same redshift, and thus (N(zf,Mf)) is largely equal for 
these objects. 

Similarly, the impact of quasar heating on low-mass, M ^ 
10 '^M0 objects has a very different dependency on the 
feedback parameter Ck than at higher mass scales. While 
(A^(0, IO'^Mq)) changes from 0.05 to 0.3 as ek goes from 0.025 
to 0.10, (A^(0,M ^ 10^3^0)) changes only from 2 to 4 for the 
same Ck values. Again this is a manifestation of feedback ef- 
fects, which are almost nonexistent in the 10 '"Mq case, but 
greatly reduce the number of outflows generated by large galax- 
ies. In fact, the true dependence of (^N(0,M ^ IO'^Mq)) on the 
feedback parameter is Ukely to be less, as eq. (33) applies the 
same correction to the luminosity function, regardless of the 
overall environment. In reality, ^(z,Lb) should be somewhat 
more suppressed in clusters than the field, leading to an even 
weaker dependence of cluster preheating on ek- 

To explore these dependencies further, in Figure 4 we plot the 
contribution to {Nsizf,Mf)) from sources above some higher 
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redshift z,, in our fiducial model. Here we see that in the higher- 
Zf, lower-mass objects, the impact of quasar halos is extremely 
weak at all times. One the other hand, in the z/ = 0, lower-mass 
objects, a substantial upturn in the contribution from sources at 
lower redshifts is seen, particularly in the IO^^Mq case. This is 
due to the ability of 5 > Scm winds from larger objects to move 
into underdense regions if given sufficient time to expand. Fi- 
nally, at higher masses, a sizeable contribution to {Nsizf,Mf)) 
is observed at all redshifts, arising from progenitor objects. 

From the values shown in Figure 3, we can directly calcu- 
late Zex(M) as well as the fraction of suppressed galaxies, as 
computed from eq. (28). These are shown in Figure 5, again 
for three representative Ck values. This plot provides a differ- 
ent perspective on the trends seen above. In our fiducial model, 
because fs(xf,M < lO'^M©) < 0.25 at all redshifts, low-mass 
galaxy formation continues largely unimpeded to this day. In 
the IO'^Mq to 1O"M0 mass range, however, a steep rise in 
Zex take place as feedback becomes more effective. Finally, at 
masses ^ lO'^-^M© , Zex values are largely independent of mass, 
and strong feedback occurs in all cases. These trends persist in 
the Ck = 0.025 and Ck = 0.10 models although the overall val- 
ues shift somewhat. In general, raising the feedback efficiency 
causes the transition region in Zex to move to smaller masses, 
while pushing the high-mass Zex plateau to higher redshifts. 

5.2. Entropy and the Quasar Luminosity Function 

The trends seen in Figure 5 have a direct impact on the B- 
band quasar luminosity function, as shown in Figure 6. From 
this plot it is iimnediately clear that the exclusion of z < Zex 
objects in our fiducial models provides an excellent fit to the 
data. Less clear is just how much of this agreement is due to 
our allowing ourselves a specific choice of Ck = 0.05. 

We would maintain that the best answer to this question is 
exactly half. There are two characteristic values that appear in 
this figure, the first of which is the characteristic luminosity, or 
mass, above which feedback effects are imprinted. This corre- 
sponds to the transition region in Figure 5, which we saw could 
be simply shifted to higher (lower) values by decreasing (in- 
creasing) Ck- Yet increases in the kinetic-energy input also have 
the effect of pushing Zex to higher redshifts. 

This corresponds to the second characteristic value in Fig- 
ure 6, namely the redshift at which feedback corrections to the 
WL03 model become important. This value is not arbitrary 
within our model, but rather is fixed when e\ is chosen to repro- 
duce the characteristic Lb- Thus the ek = 0.025 model shown 
in this figure not only fails to halt the formation of a sufficient 
number of higher-luminosity quasars at z = 0.25, but fails to 
have any impact whatsoever at z > 1 .25 . Similarly, the Ck = 0. 10 
model both over-estimates feedback effects at z = 0.25, and de- 
viates from the WL03 model too soon. 

5.3. A Heuristic Aside 

To gain a heuristic understanding of our results, it is useful to 
explore a Soltan-type argument (Soltan 1982; Yu & Tremaine 
2002), relating the accretion history of black holes to their total 
energy release over the history of the universe. We can sim- 
ply modify the argument to consider the total entropy injected 
into the IGM, rather than the total optical light emitted. To 
do so, we note from eq. (18) that Mb,heat = ekLboi%n = 
ekEradA/i/iC^fdyn = CkSrad AMbhC^. Here, the radiative efficiency 
is Crad = L\xa/{Mc^) « 0.1 from comparing estimates of the lo- 
cal black hole mass density with the quasar luminosity function 



from the 2dF redshift survey (Yu & Tremaine 2002). Thus, the 
total mass density of gas heated above San is directly propor- 
tional to the mass density in black holes, if we ignore the time 
delay due to the finite bubble propagation speed. Specifically, 
we obtain 

PbiS > 5crit) « 10Vh5l00,cri,(z)"kl+Z)"V^' (0^3 ) (Iy) . 

■(35) 

The comoving mass density in black holes as a function of red- 
shift is 

Pbbiz)= dz' dLB^nLB,z')-^. (36) 

Jz Jo EradC^ dz' 

We can use eqs. (35) and (36) to obtain 0^(5 > 5ciit,z), with- 
out recourse to Press-Schechter theory. We can also approxi- 
mately understand the exclusion redshift Zex from the fact that 
the highest density peaks are always the first to be polluted to 
high entropy, just as they are also always the first to be pol- 
luted to high metallicity. The exclusion redshift can therefore 
be approximated by: 

fiS > 5crit) = = « erfc . (37) 

Pb \V2aiM,Zex)J 

From this we can understand the evolution of the mass suppres- 
sion scale. Often in models of structure formation, the sup- 
pression mass (for instance, the Jeans mass) is self-regulating 
and therefore traces the non-Unear mass scale (e.g., Chiu & Os- 
triker 2000). This is manifestly not the case for preheating. For 
instance, if preheating affected only ^ 2(t perturbations, then 
fiS > 5ciit) ~ 0.05, independent of redshift. Instead, from the 
sharp rise in the black hole mass density with time, we know 
that f(S > 5crit) increases with time. Thus, at late times, pre- 
heating affects progressively smaller perturbations: from ~ 6a 
perturbations at z ~ 5 (Mjuppress ~ IO'^-^Mq) to ^ la perturba- 
tions at z ^ 1 (A/supress ^ 10^^ Mq). Once preheating affects the 
bulk of halos, the supply of cold gas is cut off and the comov- 
ing luminosity density in quasars and star formation abruptly 
drops. 

5.4. Entropy and Galaxy Formation 

In Figure 7 we compare our feedback models with what has 
emerged as one of the most widely-used quantifiers of galaxy 
evolution, the comoving space density of the star formation 
rate, p*. In this plot, the barrage of multi-wavelength observa- 
tions illustrates that rises by an order of magnitude between 
z = and z = 1 , and evolves much more slowly at higher red- 
shifts. For reference, we also include the analytical model given 
by eq. (2) of Hernquist and Springel (2003), with a normaliza- 
tion of p*(0) = 0.04Mq yr-i Mpc'^ 

While it is still unclear whether the evolution of reaches 
a peak around z « 1.5 (e.g. Madau et al. 1996) or stays flat to 
much higher redshifts, the primary uncertainties in this compar- 
ison are theoretical rather than observational. A rough estimate 
of the number of stars resulting from starbursts can be taken 
by simply evaluating the time derivative of the total mass con- 
tained in galaxies, scaled to the average number of stars formed 
in each burst. This gives 

il,h dz 1"°° dn r n 

Mz) = A FT / dM — — 1 -fsiz,M)\ , (38) 

^0 dt Ju^,,) dMdz 
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where is a "star formation efficiency" that parameterizes the 
fraction of gas in a given starburst that is converted to stars. 
Taking to be a free parameter that we match to observa- 
tions at high redshifts yields the soUd line shown in this figure. 
Obviously the fall-off at lower redshifts is much more severe 
than observed, yet the reason for the discrepancy is clear, as the 
majority of observed low-redshift star formation occurs not as 
bursts, but rather though quiescent star formation in disk galax- 
ies (e.g. Tan, Silk, & Balland 1999). 

This sort of star formation could simply be modeled by a 
term that depends on the total (integral) mass in galaxies 

P*,quiescent(z) = /*,?— (2)7^ / dz' X 
Vc f2o J, 

JM^oiz') dMdz' 

where the dynamical time (rnjvc) = 30fdyn as given by eq. (1) 
and is now a parameter that determines the average rate 
of quiescent star formation in disks. Adding this contribution 
(with = 0.003) to eq. (38) results in the dashed line shown in 
Figure 7. While this provides a reasonable fit to observations, 
the necessary introduction of a second free parameter makes 
this comparison a bit unsatisfying. We can conclude only that 
/9+(z) measurements are consistent with strong-feedback mod- 
els, although they do not provide particularly strong constraints. 
We note in passing, however, that our fiducial Ck = 0.05 model is 
consistent with observations of Lyman-break galaxies at z ^ 3, 
as only halos of M ;> IO^^Mq have Zex > 3, while Lyman-break 
galaxies are thought to reside in halos with masses ^ IO'^Mq, 
or even less (e.g. Scannapieco & Thacker 2003). 

A more rigorous confrontation with observations arises from 
restricting our attention to the subset of galaxies in which 
quasar feedback has the strongest impact, namely those with 
masses ^ IQ^^Mq . Thus in Figure 8 we plot the number density 
of the most luminous galaxies at three representative redshifts. 
Here we restrict our attention to K, measurements as taken by 
Pozzetti et al. (2003), as such infrared measurements are most 
sensitive to the overall galaxy stellar masses, rather than a com- 
bination of stellar masses and formation histories (e.g. Madau, 
Pozzetti, & Dickinson 1998). 

To convert our mass function to /iT-band magnitudes in an 
approximate and simple manner we make use of the Tully- 
Fisher relation observed by Pierini & Tuffs (1999), which re- 
lates the /T'-band absolute magnitude and circular velocities of 
disk galaxies as 

Mk' =-9.66v,-l.4l. (40) 

For purposes of this comparison, we associate Vc in this equa- 
tion with the underlying halo circular velocity, making no at- 
tempt to correct for galaxy morphology or differences between 
K' and bands. This allows us to compute the luminosity 
function as 

, dvr dM dn r , n 

*-^^^ = i ^^^^^t^-^^^^'^^]' '''' 

where M is a function of z' that is fixed by eqs. (40) and 
(2). This is shown as the solid (ek = 0.05) and short-dashed 
(ek = 0.025; e = 0.10) lines in Figure 8. 

Here we see that our simple model provides an excellent fit 
to the data, so much so that some degree of chance must be 



involved, given the uncertainties. This agreement is even more 
striking given the large inconsistencies between the data and the 
more detailed galaxy-formation models compiled from Menci 
et al. (2002), Cole et al. (2000), and Kauffmann et al. (1999) 
in this figure, particularly for the brighter objects. Furthermore, 
this sample is by no means the only one that illustrates a dis- 
crepancy between observations and modeling. Rather, a large 
number of independent surveys are consistently pointing to a 
lack of galaxy evolution below z ^ 1 .5. (e.^. Cimatti et al. 2002; 
Chen et al. 2003a; Glazebrook et al. 2003; Somerville et al. 
2003), which appears to present a "paradox" for current hierar- 
chical models {e.g. Chen et al. 2003b). Given this well-known 
discrepancy, is it not perhaps the case that we obtain agreement 
by overlooking an important piece of galaxy-formation physics 
that when accounted for would completely change our model? 

While a key issue in galaxy formation has indeed been ne- 
glected, we would suggest that this points not to a shortcom- 
ing in our approach, but rather to its greatest strength. Despite 
their diversity, the literature models compiled in this plot share 
the same core mechanism for arriving at the characteristic lu- 
minosity (L*), above which the number densities of galaxies 
falls off precipitously. Typically this is associated with a post- 
viralization cooling condition (Rees & Ostriker 1977; Dekel & 
Silk 1986), which arises when adopting a picture in which gas 
is shocked to the virial temperature at the redshift of halo col- 
lapse, and then must have time to cool into a disk before star- 
formation can occur. This can be implemented by enforcing 
some variation on the condition 

/?cool(A/,Zcol) > ^obs - W, (42) 

where fobs is time corresponding to the redshift of observa- 
tion, and tcoo\{M ,Zco\) is the cooling time as calculated from 
eq. (14) assuming a density of 180 times the mean density 
at Zcoi and a temperature equal to the viralization temperature 
Jvir = 35/:[v,/(kms-i)]2. 

At the lowest redshifts, this condition provides a similar and 
equally good fit to the data as our modeling, as illustrated in 
Figure 9. In the upper left panel of this plot, we see that / can 
easily be chosen to exclude star-formation in a similar subset 
of halos as quasar feedback. Yet once / is fixed, the imposition 
of a post-virialization cooling condition has a very different im- 
pact on higher-redshift halos. 

Although the (1+z)^ increase in the mean density shortens 
higher-redshift cooling times, this effect is somewhat counter- 
acted by the strong f cx ( 1 + z)-'''^ dependence on time with red- 
shift. This means that in higher-redshift objects, eq. (42) con- 
tinues to enforce a substantial redshift delay (Az) between the 
halo collapse and star formation redshifts. Thus, the reduction 
in objects with luminosities ~ at z = is achieved as a result 
of a strong redshift evolution in their number densities. This 
generic feature, while lessened somewhat by more complicated 
modeling, is the primary source of the discrepancy between the 
observations and standard models. 

Furthermore, the imposition of a post-virialization cool- 
ing criterion fails quantitatively in comparison with numerical 
simulations. If no feedback or preheating is invoked, high- 
resolution simulations predict that ~ 30-40% of the gas in 
clusters should cool (Suginohara & Ostriker 1998; Lewis et al. 
2000; Dave et al. 2001). This is because while fcooi > tu at the 
virial radius, ?cooi ^ in the dense inner regions of clusters, 
resulting in catastrophic cooUng. By contrast, K-band observa- 
tions indicate a cold gas fraction of 5- 10% in most groups 
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and clusters, independent of virial temperature. A lucid recent 
discussion of this overcooling problem is given in Balough et 
al. (2001). 

On the other hand, our proposed feedback criterion differs 
from the customary one in that it is evaluated not at the virial 
radius but at the dense centers of halos (which are the most 
highly biased and therefore most likely to have been preheated). 
Quasar feedback excludes a fixed fraction of galaxies within 
halos of a given mass and formation redshift, regardless of the 
redshift at which they are observed. Thus the solid and dotted 
lines in Figure 9 remain fixed at all redshifts, evolving only to 
enforce Zcoi > Zobs at the small-mass end. This results in the 
observed lack of evolution seen in Figure 8, which naturally re- 
produces the observations using only a simple transformation 
between v, and Mk, ■ 

Our point is not that a lack of z ^ 1.5 evolution can not be ac- 
commodated in standard galaxy-formation scenarios, but rather 
that this behavior can be reproduced only by including a series 
of detailed corrections. In galaxy formation models regulated 
by quasar feedback, this lack of low-redshift galaxy evolution 
emerges naturally, imposed by the redshift at which preheating 
becomes prevalent, which in turn is fixed by the observed value 
ofL^. 

5.5. Entropy, Clusters and the Intergalactic Medium 

Clearly, the presence of quasar outflows has impUcations far 
beyond the formation rate of supermassive black holes and the 
regulation of star formation in the galaxies that surround them. 
In fact, as described above, the most-clearly observed imprint 
of non-gravitational heating on cluster scales is not on galax- 
ies themselves, but rather on the diffuse intracluster medium 
(Kaiser 1991). The implementation of strong quasar feedback 
as in Figure 3 naturally results in an increase in ICM entropy, 
placing it on a higher adiabat that prevents it from reaching a 
high central density during collapse, which in-tum decreases 
its X-ray luminosity {e.g. Tozzi & Norman 2001). Further- 
more for a fixed level of heating per gas particle, this effect is 
more prominent for poorer clusters, whose virial temperatures 
are comparable to this extra contribution. As a result, a T^/^ 
relation is established in hot systems and broken for colder sys- 
tems, reconciling observations with semi-analytical studies and 
numerical simulations by establishing a fixed entropy "floor" of 
SS-llOr'/^keV cml 

A further constraint on the state of ICM gas comes from 
measurements of the cosmic microwave background (CMB). 
As CMB photons are relatively low-energy, thermal motions 
in hot ionized regions are able to scatter them to higher ener- 
gies: increasing the number of energetic photons and reduc- 
ing the number of lower-energy ones. The magnitude of this 
shift is proportional to the Compton-y parameter, the tempera- 
ture of the gas convolved with the total column density along 
a line-of-sight, and is thus most severe along sight-fines that 
pass through the hot and relatively dense ICM. This so called 
Sunyaev-Zel'dovich effect (Sunyaev & Zel'dovich 1972) has 
been well studied theoretically {e.g. Aghanim et al. 1997; Re- 
fregier et al. 2000) and detected in a number of clusters {e.g. 
Grego et al. 2001; Reese et al. 2002). The possibiUty of detect- 
ing the Sunyaev-Zel'dovich effect from galactic and quasar out- 
flows has been explored by previous authors (Natarajan & Sig- 
urdsson 1999; Yamada et al. 1999; Aghanim, Bafland, & Silk 
2000; Majumdar, Nath, & Chiba 2001; Platania et al. 2002), 
and the high optical depth measured by WMAP indicates that 
it might even arise from very high redshift sources z > 10, 



where Compton cooling is extremely efficient (Oh, Cooray, & 
Kamionkowski 2003). 

Similarly, the S > Sent regions that preceded clusters, while 
more difficult to measure individually, are Ukely to have left 
a substantial imprint on the global structure of the microwave 
background. In the upper panel of Figure 10, we plot the mean 
number of 5 > Saii outflows impacting an arbitrary point in 
space, {Ns{z)), as estimated from eq. (22). Here we see that 
IGM heating is extensive for all three models under considera- 
tion, impacting > 50% of the mass before the z < 1 .5 epoch of 
cluster formation, even in the lowest-energy, Ck = 0.025 case. 

We can estimate the overaU CMB distortions from these re- 
gions as 

arc /■^"■" , dt 
y{z) = / dz'—kT{z')ne{z') 

nteC^ Jj. dz' 

= 2.4 X 10-^ 1^'°^ dz'"^^^^^^ {Ns{z')) , (43) 

with rkev(z) « iookeVcm^ 8.7 X 10-3(1 + zfS^' « 5.2 X 

10~^ inii°w^ 2 (1 +z)^, where we assume that mean heated re- 

lOOkeVcm'^ ^ ^ ' 

gions are heated to levels slightly higher than 5scrit and (5 = 4, 
such that 5(5^/3 ~ 45ciit(z). While y{z) is plotted in the cen- 
tral panel of Figure 10 over the full range of redshifts under 
consideration, note that only the present day value, ^(0), is 
accessible to observation. For all models, this value is well- 
below the observational limit on the total Compton distortion 
(Fixen 1996), but comparable to the ~ 10"^ estimates of the to- 
tal impact from gravitationally heated gas {e.g. Refregier et al. 
2000). Thus more detailed comparisons between future CMB 
measurements and outflow models are likely to be able to place 
strong constraints on the overall level of quasar feedback. 

In the lower panel of Figure 10 we turn our attention to even 
more diffuse regions of space, with overdensities of only a few. 
In order to estimate the impact of quasar winds of such circum- 
galactic gas, we modify eq. (24) as 

(A(s,cgm(z)) = — / dz' dMx (44) 

Pb,0 Jz JMrniz') 

Mb.min{z',(5[z',Zcgm(z)],Z,M} X 

b[z',M,Zcs^{z),M4o(z)]^nz,LB{z,M)], 

where now Zcgm(z) < z is chosen such that at a redshift z the re- 
gion under consideration has ceased to expand, but has not yet 
coUapsed into a bound structure, which can be achieved by re- 
quiring a "turn-around overdensity" associated with a collapse 
parameter 9 = tt in eqs. (26) and (27). Finally we take M4o(z) 
as a rough estimate of the the typical Jeans mass in the post- 
reionization IGM (Barkana & Loeb 2001). 

Surprisingly, only about 5% of the circumgalactic gas has 
been impacted by outflows at z ^ 2, but the fraction of heated 
gas increases to ~ 30-40% at z ~ 0. This may seem surpris- 
ing, since the quasar population peaks at z ~ 2, after which the 
comoving density of sources falls sharply. The reason is the fi- 
nite travel time of the bubbles, as given by equation (19). Thus, 
while the winds affect the sites of massive galaxy and quasar 
formation (which are highly clustered) early on, they only im- 
pact the low density IGM at late times. In particular, while 
quasar feedback has a large impact on the collapse of non-finear 
structures, its impact on the ^ ~ 5, z ^ 2 regions associated with 
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the Lyman-a forest is expected to be small. The feedback as- 
sumed in our modeling is therefore completely consistent with 
the observed properties of this material, although it remains that 
specialized analyses may be able to detect (or rule-out) its pres- 
ence. 

At lower redshifts, the circumgalactic gas acquires temper- 
atures that tend to be described as "warm-hot." In fact, nu- 
merical simulations have shown that by z = a large fraction 
(~ 30 - 40%) of baryons have been gravitationally heated to 
temperature of order 10^-10^ K, and that this warm-hot inter- 
galactic medium (WHIM) extends into areas with overdensities 
of only a few (e.g. Cen & Ostriker 1999; Dave et al. 2001). At 
z = Dave et al. (2001) find the temperature of this gas to be 
~ 10^-^ in (5 ~ 5 regions (see also Nath & Silk 2001). Taking 
■Sciit = 150 keV and 5 ^ 5 in the 50% of the circumgalactic gas 
impacted by quasar outflows by z ~ results in a similar value. 
Thus quasar feedback represents an energy source that is com- 
parable to gravitational heating in the WHIM and is likely to 
have a large impact on its properties. Indeed, if production of 
the WHIM via quasar outflows is important, then the majority 
of the IGM will not be graviationally shocked at late times, but 
instead evolve isentropically. While observational constraints 
on this material are relatively weak at the moment, as in the 
CMB case, future studies are likely to provide a sensitive probe 
into heating of this mutiphase gas (e.g. ShuU, Tumlinson, & 
Giroux 2003) and the impact of quasar feedback. 

5.6. Relation to Baryonic Stripping 

Finally, we relate the entropy-driven process of quasar feed- 
back to the effects of starburst-driven winds. While previ- 
ous work (Scannapieco, Ferrara, & Broadhurst 2000; Thacker, 
Scannapieco & Davis 2002) has shown that dwarf outflows 
are relatively inefficient at heating gas to 5 > 5scrit levels, they 
are nevertheless able to impede the formation of galaxies with 
masses ^ IO^Mq, through the transfer of momentum. In this 
case an SN driven shock accelerates the baryons in a neighbor- 
ing perturbation to above their escape velocity, stripping them 
entirely from their associated dark matter and halting their col- 
lapse. Setting the mass of the perturbation times its escape ve- 
locity equal to the impinging momentum gives a rough estimate 
as to when such "baryonic stripping" wiU take place: 



(45) 



where Mb.p, and are the baryonic mass and radius of the per- 
turbation, while Vp = ^GMp/(2Rp), with Mp the total mass of 
the perturbation. Relating and R, to E(,o through eq. (10) 
and solving for the (physical) radius within which baryons wfll 
be stripped yields 



/?strip = 35Mpc£60 ^m;^/^ (l+z) 



(46) 



where 5s and 5p are the average overdensity of the expanding 
shock and the perturbation respectively, and M12 is the mass 
of the perturbation in units of lO'^M©. We can therefore in- 
clude baryonic stripping in our modeling of quasar feedback 
by rewriting eq. (24) to account for the interactions in which 
^strip > ^ > ^heat- This gives 



Mb,stnp(z',Ss,z,M) =min{max [Mb,heat(l,z',A^), 



(47) 



f p,<5,/?3^p(z',M)] , 'fpb5sR^,(z,z',M)}, 



where we estimate Sp/S^ w 3. The resulting Zex curves are 
shown in the upper panel of Figure 1 1 for our fiducial (F = 1, 
Ck = 0.05) model. Here we see that, due to the strong clustering 
and high masses of most quasar sources, baryonic stripping has 
an impact only on late-forming low-mass halos, while S > Sait 
heating dominates for objects with masses ^ IO'^Mq. Further- 
more, such events occur well after the typical formation redshift 
of these smaU halos, and thus have a negligible impact on the 
quasar luminosity function as shown in the lower panels of this 
figure. Thus lower-redshift quasar feedback is completely dif- 
ferent in nature from high-redshift SN driven-starbursts, and is 
driven almost exclusively by 5 > Sent heating. 



6. DISCUSSION 

In this work, we have considered the impact of quasar out- 
flows on cosmological structure formation on all scales. By 
associating only a tiny fraction of a quasar's bolometric energy 
with an outflow of the form observed in radio jets or broad- 
absorption-line winds, we have shown that such objects can 
drastically transform our picture of the z 2 universe, resolving 
a number of long-standing issues. 

Quasar feedback ties together the densest objects at the cen- 
ter of galaxies with the most diffuse regions of intergalactic gas, 
impacting all intermediate structures much as SN explosions 
impact all environments within the interstellar medium. As out- 
flows heat large regions of the IGM to entropies above ^crit, they 
regulate their own formation, resulting in the observed fall-off 
in the number density of quasars below z « 2. On larger scales, 
quasar feedback halts the formation of M ^ lO^^M© galaxies, 
providing an alternative mechanism for regulating galaxy for- 
mation than the usually-imposed condition of post-virialization 
cooling. Furthermore, while is fixed in a standard picture 
only as a result of strong redshift evolution, this value appears 
in a feedback model as the nonlinear mass scale at the redshift at 
which quasar outflows heat the denser regions of gas to 5 > 5crit 
values. Thus the observed lack of z ^ 1.5 galaxy evolution, 
while arising only though detailed corrections in a standard pic- 
ture, is a general and unavoidable feature of a model regulated 
by quasar feedback. 

On intergalactic scales, feedback heats the denser regions of 
gas to values just above 5crit, providing a natural explanation 
for why the ICM of clusters of all sizes has been generally pre- 
heated to entropy values just above the minimum excess ob- 
servable. Furthermore, while the strong clustering of quasars 
precludes them from having a serious impact on the somewhat 
overdense z ^ 2 regions associated with the Lyman-a forest, 
their global impact may provide us with sensitive observational 
probes through Compton distortions in the CMB and the multi- 
phase structure of the (z = 0) WHIM. 

While in a self-similar, hierarchical model, preheating by 
progenitors would be independent of the final host halo mass 
scale, there are three major features of quasar outflows that 
break this self-similarity. The first is that energy injection is 

5/3 

a non-unear function of halo mass, E oc Mbh cx: M^^^^^, which 
causes the largest objects to be most impacted. Secondly, the 
presence of a critical entropy imposes a scale into the problem, 
as the number of supercritical outflows is self limiting, and thus 
does not greatly exceed one even for the largest objects. Fi- 
nally, high mass objects are more strongly clustered, and there- 
fore more likely to coUapse out of preheated gas than low mass 
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objects. A useful analogy is that of metallicity: the most over- 
dense regions will tend to be polluted to high metalhcity and 
high entropy first. 

The three-phase (McKee & Ostriker 1977) model so much 
revolutionized our understanding of the interstellar medium 
that it naturally lead to widespread interest in galaxy forma- 
tion though intergalactic shocking from starbursts (e.g. Ostriker 
& Cowie 1981) and even "quasar explosions" (Ikeuchi 1981). 
While we now believe structure formation is driven instead by 
dark-matter gravitational collapse, is there not a role for out- 
flows in structure regulation! 

Today, evidence for feedback effects from high-redshift star- 
burst galaxies is better described as overwhelming than mount- 
ing. Analyses of C IV Unes in quasar absorption spectra uncover 
an IGM that has been widely and inhomogenously enriched 
with stellar material (e.g. Cowie et al. 1995; Rauch, Haehnelt, 
& Steinmetz 1997) out to very high redshifts (Songaila 2001) 
and turbulently "stirred" at z ~ 3 (Rauch, Wallace, & Barlow 
2001). Large numbers of outflowing starbursts have been di- 
rectly detected, both in optical and infrared observations at 
3 z 4 (Pettini et al. 2001) and in optical observations of 
lensed galaxies at 4 ^ z ^ 6.5 (Frye, Broadhurst, & Benitez 
2002; Hu et al. 2002). Furthermore it is clear from both an- 
alytical and numerical studies that these winds have an impor- 
tant effect on the gas in growing pre-virialized density perturba- 
tions, accelerating their baryons above the escape velocity and 
thereby regulating the formation of further such dwarves. 

At the same time, evidence of similar feedback effects at 
low redshift is staring us in the face. At z ~ 2 a large num- 
ber of highly-clustered quasars are observed to be punomeUng 



the denser regions of space, providing a natural energy source 
for the substantial preheating detected in X-ray observations of 
galaxy groups and clusters at z ~ 0. These are the coUapsed 
gaseous halos that have failed to form into galaxies, and they 
look nothing like the post-virialization clouds posited in stan- 
dard scenarios. Can we really expect non-gravitational heating 
to be the primary driver of their properties and yet have played 
only a secondary role in the formation history of large galaxies? 

As surveys scour the sky, reconstructing the history of galaxy 
formation in our universe, they only serve to highlight the in- 
terdependence between these star-forming regions and the dif- 
fuse gas that surrounds them. Similarly, a large range of nu- 
merical studies have suggested that there is more to the ther- 
mal and accretion history of the IGM than we had first imag- 
ined. By complementing galaxy surveys with detailed mea- 
surements of Compton distortions in the CMB, probes into the 
structure of the low-redshift WHIM, and perhaps even special- 
ized analyses of the Lyman-a forest, this history is accessible to 
us. Far from an abstract conjecture, quasar-regulated structure 
formation provides a clear and testable alternative to current 
low-redshift modeling, a postulated missing link that may soon 
be unearthed by multi-wavelength excavations into our cosmic 
past. 
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Fig. 1. — Evolution of the B-band quasar luminosity function. Here the data points are taken from Fan (2001, triangles), Pei (1995, circles) as compiled from 
Hartwick & Schade (1990) and Warren, Hewett, & Osmer (1994), and Barger et al. (2003, squares) as converted from the X-ray in WL03. In all panels the solid 
lines are the simple WL03 model, with F = 1.0. 




4 5 6 7 8 



log(T) (K) 




z 



Fig. 2. — Top: The cooling time at constant entropy as a function of temperature, for Sioo = 1- Since ?cool tx S^^q, changing 5ioo only shifts the curves up or down 
the vertical axis. The cooling time has a deep minimum at T^\n = 2.3 X lO^K. If the cooling time exceeds the Hubble time at r^j„, it will never be able to cool. 
Bottom: The critical entropy Scrit for which ?coo/ (^Inin) > 'h, as a function of redshift. 
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Fig. 5. — "Exclusion redshifts" for various models of quasar heating. In each panel the solid lines show Zcx such that (A'^5(zcx, A/)) = 1, while the short-dashed 
lines correspond (from top to bottom) to z values such that fs(z,M) = 1 -exp (A'i(z,M))j = 0.25, 0.5, 0.75, and 0.9. The panels are labeled by their assumed ej^ 
values, and in all cases F = I. 
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Fig. 6. — Evolution of the B-band quasar luminosity function. In all panels the long-dashed lines are the WL03 model with F = 1.0, the solid lines are our fiducial 
F = 1.0, eic = 0.05 model, and the upper and lower short-dashed lines are models in which F = 1.0 and ek is 0.025 and 0. 10 respectively. Points are as in Figure I. 
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Fig. 7. — Observed and predicted star formation rate densities. Here the solid and sliort-dashed lines represent the = 0.1, /* ^ = 0.0 and = 0.05, = 0.003 
models, respectively, the dotted line is the fit by Hemquist & Springel (2003), and the points are taken from a wide range of optical, fai-infrared, and 1.4 GHz 
measurements, as compiled and connected for reddening by Hopkins et al. (2001), and converted to our assumed cosmology. In particular, the symbols represent 
measurements by Haarsma et al. (2000) (six-pointed stars), Hopkins, Conolly, & Szalay (2000) (solid hexagons), Sullivan et al. (2000) (five-pointed stars), Serjeant 
et al. (2002) (open pentagons), Steidel et al. (1999) (filled pentagons), Yan et al. (1999) (open squares), Treyer et al. (1998) (filled squares), Tresse & Maddox (1998) 
(open triangles), Connolly et al. (1997) (filled triangles), Rowan-Robinson et al. (1997) (open circles), Gallego et al. (1995) (filled circles), and Condon (1989) 
(crosses). 




Fig. 8. — The rest-frame band luminosity function of luminous galaxies. In each panel the solid lines represent our fiducial F = 1, ei^ = 0.05 model, while 
the £]; = 0.025 and e^^ = 0. 10 models ai'e given by the upper and lower short-dashed lines respectively. The points are the Ks band observations from Pozzetti et al. 
(2003), and the dotted, long-dashed, and dot-dashed lines represent the Menci et al. (2002), Cole et al. (2000), and Kauffmann et al. (1999) models, respectively, as 
compiled by Pozzetti et al. (2003). Quasar feedback naturally turns off massive galaxy formation at low redshifts, reproducing the observed lack off evolution in the 
number density of luminous galaxies below z < 1.5. 
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Fig. 9. — Comparison between galaxy regulation by post-virialization cooling and quasar feedback. In each panel the solid line shows the last redshift of formation 
as defined by Zc.x in our fiducial model, while the short-dashed and long-dashed lines impose post-virialization criteria according to eq. (42), with / = 1 and / = 1/2, 
respectively. Finally the dotted lines correspond (from top to bottom) to z values such that fs(z,M) = 0.25, 0.5, 0.75, and 0.9, as in Figure 5. 




Fig. 10. — Global impact of quasar outflows on the intergalactic medium. Top: The mean number of 5 > Scrii outflows impacting a random point in space 
as computed from eq. (20). Here the solid line is given by our fiducial ei^ = 0.05 model, while the upper and lower dashed fines con'espond to the ej^ = 0.10 and 
e]j = 0.025 models respectively. Center: Compton-y parameter as computed from eq. (43). Solid and dashed lines are as in the top panel, while the dotted line is 
the observational limit of y < 1.5 X 10"^ (Fixen et al. 1996). Bottom: Mass fraction of heated gas in a typical region of circumgalactic gas, associated with the Lya 
forest at z > 2 and the WHIM at low redshifts, as computed from eq. (45). Lines are as in the top panel. 
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Fig. 1 1. — Comparison between feedback through S > S^it heating and the baryonic stripping of objects, which dominates in high-redshift starbursts. Top: 
As in Figure 5, the solid line shows Zex such that {Ns(zex,M)) = 1, while the short-dashed lines correspond (from top to bottom) to z values such that fg(z,M) = 
1 — exp 1^— {Ns{z,M))^ = 0.25, 0.5, 0.75, and 0.9. While quasars are able to strip the gas out of small perturbations that are not heated above 5crit, this occurs relatively 
late, long after the majority of objects in this mass range have formed. Bottom: Quasar luminosity function. In all panels the solid line is our fiducial model, but 
including the additional feedback from baryonic stripping, while the long-dashed lines are the WL03 model. Points are in Figure 1 . Stripping occurs too late to have 
any impact on these results. 



